Interpretation of the Hubble diagram in a non-homogeneous universe 
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In the standard cosmological framework, the Hubble diagram is interpreted by assuming that the 
light emitted by standard candles propagates in a spatially homogeneous and isotropic spacetime. 
However, the light from "point sources" ~ such as supernovae - probes the universe on scales where 
the homogeneity principle is no longer valid. Inhomogeneities are expected to induce a bias and a 
dispersion of the Hubble diagram. This is investigated by considering a Swiss-cheese cosmological 
model, which (1) is an exact solution of the Einstein field equations, (2) is strongly inhomogeneous on 
small scales, but (3) has the same expansion history as a strictly homogeneous and isotropic universe. 
By simulating Hubble diagrams in such models, we quantify the influence of inhomogeneities on 
the measurement of the cosmological parameters. Though significant in general, the effects reduce 
drastically for a universe dominated by the cosmological constant. 

PACS numbers: 98.80 ± k, 04.20 ± q, 42.15 ± i. 



I. INTRODUCTION 

The standard physical model of cosmology relies on a 
solution of general relativity describing a spatially homo- 
geneous and isotropic spacetime, known as the Friedmann- 
Lemaitre (FL) solution (see e.g. Ref. P]). It is assumed 
to describe the geometry of our universe smoothed on 
large scales. Besides, the use of the perturbation theory 
allow^s to understand the properties of the large scale 
structure, as well as their growth from initial conditions 
set by inflation and constrained by the observation of the 
cosmic microwave background. 

While this simple solution of the Einstein field equa- 
tions, together with the perturbation theory, provide a 
description of the universe in agreement with all exist- 
ing data, it raises many questions on the reason why it 
actually gives such a good description. In particular, it 
involves a smoothing scale which is not included in the 
model itself ^ . This opened a lively debate on the fitting 
problem [3] - i.e. what is the best-fit FL model to the 
lumpy universe? - and on backreaction - i.e. the fact 
that local inhomogeneities may affect the cosmological 
dynamics. The amplitude of backreaction is still actively 
debated [IHS] (see Ref. [7] for a critical review). 

Regardless of backreaction, the cosmological model as- 
sumes that the distribution of matter is continuous (i.e. it 
assumes that the fluid approximation holds on the scales 
of interest) both at the background and perturbation 
levels. Indeed numerical simulations fill part of this gap 
by dealing with A^-body gravitational systems in an ex- 
panding space. The fact that matter is not continuously 
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distributed can however imprint some observations, in 
particular regarding the propagation of light with narrow 
beams, as discussed in details in Ref. [8]. It was argued 
that such beams, as e.g. for supernovae observations, 
probe the spacetime structure on scales much smaller 
than those accessible in numerical simulations. The im- 
portance of quantifying the effects of inhomogeneities on 
light propagation was first pointed out by Zel'dovich [3]. 
Arguing that photons should mostly propagate in vacuum, 
he designed an "empty beam" approximation, generalized 
later by Dyer and Roeder as the "partially filled beam" 
approach [TU] . More generally, the early work of Ref. fS] 
stimulated many studies on this issue. [TTH25] . 

The propagation of light in an inhomogeneous universe 
gives rise to both distortion and magnification induced by 
gravitational lensing. While most images are demagnified, 
because most lines of sight probe under-dense regions, 
some are amplified because of strong lensing. Lensing 
can thus discriminate between a diffuse, smooth compo- 
nent, and the one of a gas of macroscopic, massive objects 
(this property has been used to probe the nature of dark 
matter on galactic scales [Ml. Therefore, it is expected 
that lensing shall induce a dispersion of the luminosities 
of the sources, and thus an extra-scatter in the Hubble 
diagram [57]. Indeed, such an effect does also appear at 
the perturbation level - i.e. with light propagating in 
a perturbed FL spacetime - and it was investigated in 
Refs. [3S1 [211 |S5] . The dispersion due to the large-scale 
structure becomes comparable to the intrinsic dispersion 
for redshifts z > 1 [SU] but this dispersion can actually 
be corrected [311 - 136] . Nevertheless, a considerable frac- 
tion of the lensing dispersion arises from subarcminute 
scales, which are not probed by shear maps smoothed 
on arcminute scales [37| . The typical angular size of the 
light beam associated with a SN is typically of order 
10"'' arcsec (e.g. for a source of physical size ^ 1 AU at 
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redshift z ~ 1), while the typical observational aperture is 
of order 1 arcsec. This is smaller than the mean distance 
between any massive objects. 

One can estimate [38] that a gas composed of parti- 
cles of mass M can be considered diffuse on the scale 
of the beam of an observed source of size if M < 
2 X lO~^^M0/i^ ('^s/l AU)^. In the extreme case for which 
matter is composed only of macroscopic pointlike objects, 
then most high-redshift SNela would appear fainter than 
in a universe with the same density distributed smoothly, 
with some very rare events of magnified SNela 38-40 . 
This explicits the connection between the Hubble diagram 
and the fluid approximation which underpins its standard 
interpretation. 

The fluid approximation was first tackled in a very 
innovative work of Lindquist and Wheeler [32] , using a 
Schwarzschild cell method modeling an expanding uni- 
verse with spherical spatial sections. For simplicity, they 
used a regular lattice which allows a very limited set 
of possibilities, corresponding to the most homogeneous 
topologies of the 3-sphere [43]. It has recently been revis- 
ited in Refs. [33] and in Refs. [3S] for Euclidean spatial 
sections. They both constructed the associated Hubble 
diagrams, but their spacetimes are only approximate so- 
lutions of the Einstein field equations. An attempt to 
describe filaments and voids was also proposed in Ref. [3^ . 

These approaches are conceptually different from the 
solution we adopt in the present article. We consider an 
exact solution of the Einstein field equations with strong 
density fluctuations, but which keeps a well-defined FL 
averaged behaviour. Such conditions are satisfied by the 
Swiss-cheese model [33- one starts with a spatially ho- 
mogeneous and isotropic FL geometry, and then cuts out 
spherical vacuoles in which individual masses are embed- 
ded. Thus, the masses are contained in vacua within 
a spatially homogeneous fluid-filled cosmos (see bottom 
panel of Fig. [2]). By construction, this exact solution 
is free from any backreaction: its cosmic dynamics is 
identical to the one of the underlying FL spacetime. 

From the kinematical point of view, Swiss-cheese models 
allow to go further than perturbation theory, because not 
only the density of matter exhibits finite fluctuations, 
but also the metric itself. Hence, light propagation is 
expected to be very different in a Swiss-cheese universe 
compared to its underlying FL model. Moreover, the 
inhomogeneities of a Swiss-cheese are introduced in a 
way that addresses the so-called "Ricci-Weyl problem" . 
Indeed, the standard FL geometry is characterized by a 
vanishing Weyl tensor and a non-zero Ricci tensor, while 
in reality light mostly travels in vacuum, where conversely 
the Ricci tensor vanishes - apart from the contribution 
of A, which does not bend light - and the Weyl tensor is 
non-zero (see Fig.[T]). A Swiss-cheese model is closer to 
the latter situation, because the Ricci tensor is zero inside 
the holes (see Fig. [2]) . It is therefore hoped to capture 
the relevant optical properties of the universe. 

In fact, neither a Friedmann-Lemaitre model nor a 
Swiss-cheese model can be considered a realistic descrip- 



tion of the universe. They share the property of being 
exact solutions of the Einstein equations which satisfy 
the Copernican principle, either strictly or statistically. 
Swiss-cheese models can be characterized by an extra- 
cosmological parameter describing the smoothness of their 
distribution of matter. Thus, a FL spacetime is nothing 
but a perfectly smooth Swiss-cheese. It is legitimate to 
investigate to which extent observations can constrain 
the smoothness cosmological parameter, and therefore to 
quantify how close to a FL model the actual universe is. 
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FIG. 1: The standard interpretation of SNe data assumes 
that light propagates in purely homogeneous and isotropic 
space (top). However, thin light beams are expected to probe 
the inhomogeneous nature of the actual universe (bottom) 
down to a scale where the continous limit is no longer valid. 

The propagation of light in a Swiss-cheese universe 
was flrst investigated by Kantowski [48 , and later by 
Dyer and Roeder [32] ■ Both concluded that the effect, on 
the Hubble diagram, of introducing "clumps" of matter 
was to lower the apparent deceleration parameter. The 
issue was revived within the backreaction and averaging 
debates, and the Swiss-cheese models have been extended 
to allow for more generic distributions of matter inside 
the holes - instead of just concentrating it at the center - 
where spacetime geometry is described by the Lemaitre- 
Tolman-Bondi (LTB) solution. The optical properties 
of such models has been extensively studied (see Refs. 
[SOHSH]) to finally conclude that the average luminosity- 
redshift relation remains unchanged with respect to the 
purely homogeneous case, contrary to the early results of 
Refs. [351I3S1- 

In general, the relevance of "LTB holes" in Swiss-cheese 
models is justified by the fact that they allow to reproduce 
the actual large-scale structure of the universe (with voids 
and walls). However, though inhomogeneous, the distri- 
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FIG. 2: Swiss-cheese models (bottom) allow to model inhomo- 
geneities beyond the continous limit, while keeping the same 
dynamics and average properties as the FL model (top). 



button of matter in this class of models remains continous 
at all scales. On the contrary, the old-fashioned approach 
with "clumps" of matter inside the holes breaks the contin- 
uous limit. Hence, it seems more relevant for describing 
the small-scale structure probed by thin light beams. 

In this article, we revisit and update the studies of 
Refs. [m |49] within the paradigm of modern cosmology. 
For that purpose, we first provide a comprehensive study 
of light propagation in the same class of Swiss-cheese 
models, including the cosmological constant. By gen- 
erating mock Hubble diagrams, we then show that the 
inhomogeneities induce a significant bias in the apparent 
luminosity-redshift relation, which affects the determina- 
tion of the cosmological parameters. As we shall see, the 
effect increases with the fraction of clustered matter but 
decreases with A. For a universe apparently dominated 
by dark energy, the difference turns out to be small. 

The article is organized as follows. Section [H] describes 
the construction and mathematical properties of the Swiss- 
cheese model. In Section |III[ we summarize the laws of 
light propagation, and introduce a new tool to deal with a 
patchwork of spacetimes, based on matrix multiplications. 
In Section |IV[ we apply the laws introduced in Sec. |III| to 
Swiss-cheese models and solve the associated equations. 
The results allow to investigate the effect of one hole (Sec- 
tion |v|) and of many holes (Section |VI[ ) on cosmological 
observables, namely the redshift and the luminosity dis- 
tance. Finally, the consequences on the determination of 



II. DESCRIPTION OF THE SWISS-CHEESE 
COSMOLOGICAL MODEL 

The construction of Swiss-cheese models is based on 
the Einstein-Straus method [17] for embedding a point- 
mass within a homogeneous spacetime (the "cheese" ) . It 
consists in cutting off a spherical domain of the cheese 
and concentrating the matter it contained at the center of 
the hole. This section presents the spacetime geometries 
inside and outside a hole (Subsec II A), and how they are 
glued together (Subsec. II B). 



A. Spacetime patches 

1. The "cheese" - Friedmann-Lemaitre geometry 

Outside the hole, the geometry is described by the 
standard Friedmann-Lemaitre (FL) metric 

ds^ = -dT' + a\T) [dx' + fUx) df^'] , (2.1) 

where a is the scale factor and T is the cosmic time. The 
function fxix) depends on the sign of K and thus of the 
spatial geometry (spherical, Euclidean or hyperbolic), 
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respectively for K > K = Q ox K < Q. The Einstein 
field equations imply that the scale factor a{T) satisfies 
the Friedmann equation 
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and where p = po(ao/a)'^ is the energy density of a pres- 
sureless fluid. A subscript indicates that the quantity 
is evaluated today. It is convenient to introduce the 
cosmological parameters 
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in terms of which the Friedmann equation takes the form 

|-)'=an(^)V^^.(^)Vf^A. (2.5) 



2. The "hole" - Kottler geometry 

Inside the hole, the geometry is described by the exten- 
sion of the Schwarzschild metric to the case of a non-zero 
cosmological constant, known as the Kottler solution 
[5^ [50] (see e.g. Ref. [51] for a review). In spherical 
coordinates {r,6,(p), it reads 



the cosmological parameters is presented in Section VII 



ds^ = -A{r) dt^ + A-\r) dr^ + r^dft^, 
with A{r) = 



r 



(2.6) 
(2.7) 
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and where rs = 2 GM is the Schwarzschild radius asso- 
ciated with the mass M at the center of the hole. It is 
easy to check that the above metric describes a static 
spacetime. The corresponding Kilhng vector = Sq has 
norm g^uS.^S,'^ = A(r) and is therefore timehke as long as 
A> Q. Hence, there are two cases: 

1. If %GMf^ > 1, then A{r) < for all r > 0, 
so that is spacelike. In this case, the Kottler 
spacetime contains no static region but it is spatially 
homogeneous. 

2. If 9{GM fA < 1, then A{r) > for r between 
and Tc > Tb which are the two positive roots of the 
polynomial rA{r), and correspond respectively to 
the black hole and cosmological horizons. We have 
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(2.10) 
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with cos ij — 3 GM\/A, so that 

3 1 3 

fs < ''b < < -7= < f c < -T=- 
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In the region rb < r < rc, the Kottler spacetime is 
static. Note also that r = and r = are Killing 
horizons, since ^ vanishes on these hypersurfaces. 

In practice, we use the Kottler solution to describe the 
vicinity of a gravitationally bound object, such as a galaxy, 
or a cluster of galaxies. In this c ontext , we have typically 
9(GM)2A < 10"^** (see Subsec. |V A[), so we are in the 



second case. Moreover, this solution only describes the 
exterior region of the central object, it is thus valid only 
for r > Tphys, where ^phys is the physical size of the object. 
For the cases we are interested in, ^phys 3> ^b, so that 
there is actually no black-hole horizon. 



B. Junction conditions 

Any spacetime obtained by gluing together two different 
geometries, via a hypersurface E, is well-defined if - and 
only if - it satisfies the Israel junction conditions [62l l63]: 
both geometries must induce (a) the same 3-metric, and 
(b) the same extrinsic curvature on E. 

The junction hypersurface E is the worldsheet of a 
comoving 2-sphere, as imposed by the symmetry of the 
problem. Hence, it is defined hy x = Xh = est in FL 
coordinates, and by r = ni(i) in Kottler coordinates. 
Both points of view are depicted on Fig. |3] 

In the FL region, the normal vector to the hypersurface 
is given by n^ff^'^ — S^/a. The 3-metric and the extrinsic 
curvature induced by the FL geometry are respectively 
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FIG. 3: The junction hypersurface as seen from the FL point 
of view with equation X = Xh (left); and from the Kottler 
point of view with equation r — rh{t) (right). 



where (x") — (T, 6, tp) are natural intrinsic coordinates 
for E. We stress carefully that, in the following and as 
long as there is no ambiguity, a dot can denote a time 
derivative with respect to T or t, so that a — da/dT and 
Th = dr^i/dt, while a prime can denote a derivative with 
respect to x or r, so that — df^/dx and A' = dA/dr. 
The 3-metric induced on E by the Kottler geometry is 



ds2 = -K^{t)df + rl{t)dn^, 



where 



lA^r^{t)]-rl{t) 



Therefore, the first junction condition implies 

rh(t) = a(T)xh, 
dT 

d^^'^W' 



(2.13) 

(2.14) 

(2.15) 
(2.16) 



which govern the dynamics of the hole boundary, and 
relate the time coordinates of FL and Kottler regions. 

The extrinsic curvature of E induced by the Kottler 
geometry, but expressed in (x°) coordinates, reads 



K^^^dx^dx" = - 



rh + KA'(rh)/2^^2^rhA(rh)^^2 

(2.17) 



Hence, the second junction condition is satisfied only if 
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(2.18) 

It is straightforward to show that Eq. ( 2.18| ), together 
with the Friedmann equation (2.3), implies that Kottler 



and FL regions have the same cosmological constant, and 
47r 
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(2.19) 
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C. Summary 

Given a FL spacetime with pressureless matter and a 
cosmological constant, a{T) is completely determined 
from the Friedmann equation. A spherical hole of 
comoving radius Xh, which contains a constant mass 
M = 47rpa'^/|-(xh)/3 at its center, and whose geometry 
is described by the Kottler metric, can then be inserted 
anywhere. The resulting spacetime geometry is an exact 
solution of the Einstein field equations. 

By construction, the clump inside the hole does not 
backreact on the surrounding FL region. It follows that 
many such holes can be inserted, as long as they do not 
overlap. Note that if two holes do not overlap initially, 
then they will never do, despite the expansion of the 
universe, because their boundaries are comoving. 



III. PROPAGATION OF LIGHT 



A. Light rays 

The past lightcone of a given observer is a constant 
phase hypersurface w = const. Its normal vector = 
df^w (the wave four- vector) is a null vector satisfying the 
geodesic equation, and whose integral curves (light rays) 
are irrotational: 



k^k^ — 0, 



k'^V.k^ = 0, 



V[^fc^] = 0. 



(3.1) 



For an emitter and an observer with respective four- 
velocities u^j^ and u^^jjg, we define the redshift by 

, _ ^^111^^(^0111) 



-'obs"'A' 
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^'^ = da;^/d7. Hence, this vector encodes the whole 
information on the size and shape of the bundle. 

Having chosen v = at the observation event - which is 
a vertex point of the bundle - ensures that the separation 
vector field is everywhere orthogonal to the geodesies, 
fc^^^ = 0. In such conditions, the evolution of with v 



is governed by the geodesic deviation equation 

k'^k^V^Vp^'^ = R^^^^pk'^k^iP, 
where R^yaB is the Riemann tensor. 



(3.4) 



2. Sachs equation 

Consider an observer with four- velocity u^. In view 
of relating ^'^ to observable quantities, we introduce the 
Sachs basis (s^)yi6{i,2}i defined as an orthonormal basis 
of the plane orthogonal to both and fc'^. 



s^^SB^ = s'^Uf_, ^ s'^kf, = 0, (3.5) 



and parallel-transported along the geodesic bundle. 



0. 



(3.6) 



The plane spanned by (si, S2) can be considered a screen 
on which the observer projects the light beam. The 
2- vector of coordinates S^a — (.iiS^ then represents the 
relative position, on the screen, of the light spots corre- 
sponding to two neighbouring rays separated by 

The evolution of , with light propagation, is deter- 
mined by projecting the geodesic deviation equation (3.4) 
on the Sachs basis. The result is known as the Sachs 
equation [U [Ml , and reads 



(3.7) 



where v is an afhne parameter along the geodesic, so that 
k^ = dx'^/dv, and w = at the observation event. The 
wave four-vector can always be decomposed into temporal 
and spatial components. 



#d^ = l, (3.3) 



k^' = (1 -Fz)(u^ - d^') 

where d'^ denotes the spatial direction of observation. In 
Eq. (3.3), we have chosen an affine parameter adapted 



to the observer, in the sense that 27ri'o = Wobs^A'(O) = 1- 
This convention is used in all the remainder of the article. 



where TZab — Rfj.i/apk'^ k^nj^n^ is the screen-projected 
Riemann tensor, called optical tidal matrix. It is conve- 
niently decomposed into a Ricci term and a Weyl term as 
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(3.8) 



Ic^.apa^'k'^a^k^, (3.9) 



B. Light beams 

1. Geodesic deviation equation 

A light beam is a collection of light rays, that is, a 
bundle of null geodesies {a;^(w,7)}, where 7 labels the 
curves and v is the affine parameter along them. The 
relative behaviour of two neighbouring geodesies x'^(-,7) 
and x'^{-,j + d'y) is described by their separation vector 



3. Notions of distance 

Since the light beam converges at the observation event, 
we have £,^{v = 0) = 0. The linearity of the Sachs 
equation then implies the existence of a 2 x 2 matrix 
I?'^^ , called Jacobi matrix, such that 

r(«) = i?^(«)(^) (3.10) 
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From Eq. (3.7), we immediately deduce that this matrix 



satisfies the Jacobi matrix equation 
d2 



with initial conditions 
2)^b(0)=0, 



dv 



(3.12) 



We shall also use the sho rt-ha nd notation ^ = and 
T> = (V^b) so that Eq. ( |3.1l[ ) reads d'^V/dv^ V, 
with I>(0) = and i>(0) = 1. 

Since the Jacobi matrix relates the shape of a light 
beam to its "initial" aperture, it is naturally related 
to the various notions of distance used in astronomy 
and cosmology. The angular distance Da is defined by 
comparing the emission cross-sectional area d^S'sourcc of a 
source to the solid angle dfi^i^g under which it is observed. 



d'Ss. 



It is related to the Jacobi matrix by 
Da = V|detI>(wsourc 



(3.13) 



(3.14) 



where Wsourco is the affine parameter at emission. 

The luminosity distance Z^l is defined from the ratio 
between the observed flux -Fobs a-nd the intrinsic luminos- 
ity ^source of the source, so that 



(3.15) 



It is related to the angular distance by the following 
distance duality law 



Di^ = {1 + z)^Da. 



(3.16) 



Hence, the theoretical determination of the luminosity 
distance relies on the computation of the Jacobi matrix. 



C. Solving the Sachs equation piecewise 

Since we work in a Swiss-cheese universe, we have to 
compute the Jacobi matrix for a patchwork of spacetimes. 
It is tempting, in this context, to calculate the Jacobi 
matrix for each patch independently, and then try to 
reconnect them. In fact, such an operation is unnatural, 
because the very definition of "D imposes that the inital 
condition is a vertex point of the light beam. Thus, 
juxtaposing two Jacobi matrices is only possible at a 
vertex point, which is of course too restrictive for us. 

We can solve this problem by extending the Jacobi 
matrix formalism into a richer structure. This requires 



to consider the general solution of Eq. (3.7), for arbitrary 
initial conditions. Thus, we have 



(3.11) ^{v) = Civ; «init)-|.=..„.,+I>(f ; finit) 



dv 



(3.17) 



as for any linear second order differential equation, solved 
from Uinit to V. In the following, C(u;'i;init) is referred 
to as the scale matrix. It is easy to check that both 
the scale and Jacobi matrices satisfy the Jacobi matrix 



equation (3.11) but with different initial conditions: 

dV 



f (t'init;Wmit) = 0, 

whereas 

C(winit; I'init) = 1, 



dv 



dC 

dv 



(Winit;i'init) = 1, (3.18) 



(winit;Winit) = 0. (3.19) 



The most useful object for our problem turns out to be 
the 4x4 Wronski matrix constructed from C and f , 



^AU \-( C(W Winit) I>(w Winit) 

W(-y;winit)= dC / \ dx>V \ 

> — (w;Winit) -J^rl^S^init) 



dv 



(3.20) 



in terms of which the general solution (|3.17|) reads 
(v) = W(u;ui„it) 



dv 



Av 



(«init). (3.21) 



It is clear, from Eq. (3.21), that W satisfies the relation 
W{vr,v^) = W{v^;v2) ■ W(«2; «3). (3.22) 

Hence, the general solution of the Sachs equation in a 
Swiss-cheese universe can be obtained by multiplying 
Wronski matrices, according to 



W(z;s, 



; 0) = >VFL(f source; >Vk(^': 



„(1).„,(1) 



in ' out 



WFL(«;iul; «L'^) • • • y^Fi.{v)::>;Q) (3.23) 



{N). 



where Wfl a-nd Wk are the Wronski matrices computed 
respectively in the FL region and in the Kottler holes. v[ll 

and are the values of the affine parameter respectively 
at the entrance and the exit of the ith hole. 



IV. INTEGRATION OF THE GEODESIC AND 
SACHS EQUATIONS 

Consider an observer lying within a FL region, who 
receives a photon after the latter has crossed a hole. In 
this section, we determine the light path from entrance 
to observation by solving the geodesic equation, and we 
calculate the Wronski matrix for the Sachs equation. 



FIG. 4: A light ray propagates alternatively in FL and Kottler regions. The main geometrical quantities defined and used in 
Sec. |IV| are depicted on this simplified view of a single hole. 



The main geometrical quantities are summarized on 
Fig. in (j ^ is the direction of observation as defined in 
Eq. (3.3 1. The spatial sections of the FL region can 



be described either by comoving spherical coordinates 
(x, 0, (f) or, when the spatial sections are Euclidean, by 
comoving Cartesian coordinates {X, Y, Z) . 

The hole is characterized by its comoving spatial posi- 
tion X^, in terms of the FL coordinates, and its mass AI, 
or equivalently its comoving radius ph- Note that, con- 



trary to Subsec. II A it is no longer denoted Xh, in order 
to avoid confusion with the radial comoving coordinate 
of the center of the hole. 

A photon enters into the hole with wave vector fc^^, exits 
from it with wave vector A:^^^, and reaches the observer 
with wave vector fcg. We respectively denote £in, £out 
and £o the associated events. The coordinates of the 
first two can be expressed either with respect to FL, e.g. 
as (Tin,X?j) in Cartesian coordinates, or with respect to 
the hole, e.g. as (^in, Hnj ^im (ySin) in the Kottler spherical 
coordinate system. 

Our calculations go backward in time. Starting from Eq, 
we first determine fout, WFL(wout; t'obs), and second fin, 
VVK(^^in; ^^out)- The same operations can then be repeated 
starting from fin and so on. 



A. Friedmann-Lemaitre region (from £o to £out ) 

The geometry of the Friedmann-Lemaitre region is 
given by the metric (2.1 ) which can be rewritten in terms 



of the conformal time -q, defined by drj — dT/a(T), as 



ds^ = a'iv) [^drf + dx' + fUx) dn"] 



(4.1) 



1. Geodesic equation 

If one chooses the center x = of the FL spherical 
coordinate system on the worldline of the (comoving) 
observer, then the geodesic equation is easily solved as 



(4.2) 



which corresponds to a purely radial trajectory. Note 
however that for a generic origin, this is no longer true. 
The associated wave vector remains colinear to the ob- 
served one, /cq . It is only subject to a redshift induced by 
the cosmic expansion, so that 



(4.3) 



We stress that, in Eq. (|4.3|), fi = refers to components 
on d^, not on 



d,j/a. 



2. Intersection with the hole 

Once the geodesic equation has been solved and the 
position of the hole has been chosen, we can calculate 
the intersection font between the light ray and the hole 
boundary. In the particular case of a spatially Euclidean 
FL solution {K — 0), the Cartesian coordinates X^^^^ of 
£out satisfy the simple system of equations 



{K 



X, 



^out =^0 + (%-%ut)rf' 



(4.4) 



where X\^ and Xq 



are the respective Cartesian coordi- 
nates of the hole and the observer, while c?* is the spatial 
direction of observation. Although conceptually similar. 
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the determination of Sout for a FL solution with arbitrary 
spatial curvature is technically harder. 



In general, we deduce from Eq. (4.3) that the wave vec- 
tor at fout is A:^ut = (ao/aout)^^o ' where Oout = a(?7out)- 



3. Wronski matrix 

In the FL region, the Sachs basis (si,S2) is defined 
with respect to the fundamental observers, comoving with 
four- velocity u = dr- The explicit form of this basis does 
not need to be specified here. 

The Sachs equation can be solved analytically by means 
of a conformal transformation to the static metric 

dl2 ^ al [-drj^ + /i(x)dr!2] ^ ^^^dxMx". (4.5) 

Because the geometries associated with g^^ and g^i^ are 
conformal, any null geodesic for g^i, afhnely parameterized 
by V is also a null geodesic for g^^ affinely parameterized 
by -D, with a^dv — a^dv. As dv = {a^/aQ)dr], it follows 
that V = aoTj. 

For the static geometry, the optical tidal matrix reads 
TZ = —{K/aQ)l, so that the Sachs equation is simply 



dr/2 



(4.6) 



We then easily obtain the Jacobi and scale matrices: 
V = aafK'iv - Vinit)!, C = /^^ (77 - r;i„it)l. (4.7) 

To go back to the original FL spacetime, we use that 
dv = a^dr] and the fact that the screen projections of 
the separation vectors for both geometries are related by 
fl^ = tio^- The final result is 



f FL = flinit —fKiV - ?7init)l, 

CfL = /k(7? - ?7init) - T-LmitfKiV - Vinit) 

Oinit L 



(4.8) 



(4.9) 



where H = a'{jf)/a{vi) is the conformal Hubble function. 
This completely determines Wfl- 

Note that we can recover the standard expression of 
the angular distance by taking the initial condition at the 
observer. The relation (3.14) then implies 



-Da = VdetX>FL - ^KiVs. 

(1 + z) 



(4.10) 



where z = oo/a — 1 is the redshift of a photon that only 
travels through a FL region. 



B. Kottler region (from fout to £inj 

1. Initial condition at EoMt 

In the previous section, we have determined £out and 
/Cqu^ in terms of the FL coordinate system. However, in 



order to proceed inside the hole, we need to express them 
in terms of the Kottler coordinate system (t, r, 0, ip). 

A preliminary task consists in expressing £aut a-nd fc^^^ 
in terms of FL spherical coordinates, with origin at the 
center of the hole. This operation is straightforward. The 
event £out is then easi ly con verted, since (a) we are free to 
set tout = 0, (b) Eq. \2.1b\ implies rout = a('7out)Ph, and 
(c) the angular coordinates 6'out, V'out remain unchanged 
if the Kottler axes are chosen parallel to the FL ones. 

The first junction condition ensures that light is not 
deflected when it crosses the boundary E of the hole. 
Indeed, the continuity of the metric implies that the 
connection does not diverge on E. Integrating the geodesic 
equation dfc'^ = — r[^^fc"fc^dv between and v^^^ then 
shows that A:^ is continuous at fout- Therefore, we just 
need to convert its components from the FL coordinate 
system to the Kottler one. The result is 



'*'out 



- 

"-out ~ 

- 



aout 



aout 

"-out 

k'P 
"-out- 



fcout + - ^('^out)fc?u 



Vl - A(rout)Ct + ^ou 



(4.11) 

(4.12) 

(4.13) 
(4.14) 



2. Shifting to the equatorial plane 

Since the Kottler spacetime is spherically symmetric, it 
is easier to integrate the geodesic equation in the equato- 
rial plane 9 = 7r/2. In general, however, we must perform 
rotations to bring both fout and /cout into this plane. 

Starting from arbitrary initial conditions (foutj^out); 
we can shift to the equatorial plane in two steps. In the 
following, Ri{d) denotes the rotation of angle d about 
the x'-axis. The operations are depicted on Fig. [5] 

• First, bring fout to the point fout.cq on the equato- 
rial plane by the action of two successive rotations, 
R7X~Vo\it) followed by Ry{'n /2 — 6'out)- The wave 
vector after the two rotations is denoted fcout- 

• Then, bring fcout to the equatorial plane with 

(—■(/'), where is the angle between the pro- 
jection of fcout on the yz-plane and the y-axis. Note 
that such a rotation leaves fout.cq unchanged. 

It follows that, after the 3 rotations 

R = R^{-^) °Ry{l' 9out) o R.i-ifout), (4.15) 

£out and fcout are changed into £out,cq and fcout.oq which 
lie in the equatorial plane. In the following, we omit 
subscripts 'eq', keeping in mind that we will have to 
apply R^^ to recover the original system of axes. 
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"-oiil , 




FIG. 5: An arbitrary initial condition is rotated so that the geodesic lies in the equatorial plane 9 = n/2. Left: £out is brought 
(i) to (/3 = by the rotation Rz{—'fiout), and (ii) to ^ = tt/2 by the rotation Ry{n/2 ~ 6out)- The resulting event and wave vector 
are denoted 5out,cq and k'^^. Middle: k'^^ is brought to the equatorial plane by the rotation ii^; (—?/>). Right: final situation. 



3. Null geodesies in Kottler geometry 

In the Kottler region, the existence of 2 Killing vectors 
associated to statisticity and spherical symmetry implies 
the existence of 2 conserved quantities, the energy E and 
the angular momentum L of the photon. It follows that 
a null geodesic is solution of 



av 



Mr] 



L = E\ r 



.dip 
Av 



= L. 



(4.16) 

Introducing the dimensionless variable u = rg/r and the 
impact parameter b = L/E, Eqs. (4.161 imply 



du 
dip 

du 
£;2 V cb 



P{u), 



(4.17) 
(4.18) 
(4.19) 



with 



A{u) = 1 - u + £2U~'^ , P{u)=£i-u^A{u), (4.20) 

and where Si = rs/b and £2 = Arg/3. 

Our purpose is now to compute the coordinates 
itin,rin, ifiin) and the components fc(^j of the wave vec- 
tor at the entrance event given those at Sout- The 
situation is summarized on Fig. |6] 

The radius (or alternatively ujn) and time ijn at en- 
trance are determined by comparing the radial dynamics 
of the photon, governed by Eq. (4.17), to the one of the 



hole boun dary. The latter is obtained from Eqs. (2.14) 
and (2.18). By introducing Uh = rg/r^^, it reads 

duh 



rs- 



dt 



ugA(wh)Vl-^K). (4.21) 




FIG. 6: Null geodesic in the Kottler region. Depicted with 
the Kottler coordinate system, the hole grows so that the ray 
enters with rin and exits with rout > ^in. 



Eqs. (4.17) and (4.21) are then integrated"^ as iphoton(w) 
and thole (uh). The entrance radius then results from 
solving numerically the equation tphoton(uin) ~ thoio('"in), 
which also provides tin- 

The usual textbook calculation of the deflection angle 
A.ipao of a light ray in Kottler geometry yields 



2eu 1 



^ AGM 



A 62 
~3~ 



(4.22) 



at lowest order in £1 and £2. However, we cannot use this 
expression here - although it gives its typical order of 
magnitude ~ because A(^oo represents the angle between 
the asymptotic incoming and outgoing directions of a ray. 



^ The integration can be perfor med either numerically, or analyti- 
cally in the case of Eq. ||4.17[l and perturbatively for Eq. (|4.21[l . 
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whereas we must take into account the finite extension of 
the hole (see Fig. |6]). 

In general, the deflection angle Aip — (ySout ~ Vin is 



du 



du 



27r 



(4.23) 



where P{u) is the polynomial defined in Eq. (4.201, and 



1*111 is the value of u at minimal approach. The integral 
involved in Eq. (4.231 can be rewritten as 



du' 



\/P{u') V"3 - "2 



F 



U2 — U U2 — Ui 



U2 



Ui U2-U3 

(4.24) 

where Ui < U2 = < are the three (real) roots of 
P{u), and F{ip, e) denotes the elliptic function of the first 
kind 1661 



F(^,e) 



d6 



VT" 



(4.25) 



e sm 



Thus, Eq. (4.24) provides an exact expression of the de- 



flection angle A 1^9, and therefore of (/Jj, 

Once fin is determined, it is easy to obtain fc^^^ by using 
the constants of motion. The result is 



E 



A(riii) 



A{nn) inn klY 



(4.26) 
(4.27) 
(4.28) 



4- Final conditions at £in 



The last step consists in coming back to the original 
FL coordinate system. That means (a) using to 
recover the initial system of axes, and (b) converting the 
components of Sin and fc^^ in terms of the FL coordinate 
system. We have already described such operations in 



IV B 2 and 



IV B 1 respectively, except for the time 
coordinate (since we set iout = 0)- 

The easiest way to compute the cosmic time Tin at 
entrance is to use the relation r-m = a(Tln)/A'(pii). In a 
spatially Euclidean FL spacetimc {K = 0), we get 



: argsinh 



1 ~ ^Ao \aoPhJ 



(4.29) 

With this last result, we have completely determined the 
entrance event 



5. Sachs basis and optical tidal matrix 

Once the geodesic equation is completely solved, we 
are ready to integrate the Sachs equation in the Kottler 



region, that is, to determine the Wronski matrix Wk- 
Such a task requires first to define the Sachs basis (si, S2) 
with respect to which Wk will be calculated. 

The four-velocity u is chosen to be the one of a radially 
free-falling observer. 



1 



Mr) 



dt + ^l-A{r)dr. 



(4.30) 



This choice ensures the continuity of u through the hole 
frontier, where u = Dt- The wave 4- vector k is imposed 
by the null geodesic equations, and reads 



k^ -^dt±E\ll 
A{r) 



^ J.2 



(4.31) 



where the ± sign depends on whether the photon ap- 
proaches (— ) or recedes (-I-) from the center of the hole. 

By definition, the screen vectors si, S2 form an or- 
thonormal basis of the plane orthogonal to both u and k. 
Here, since the trajectory occurs in the equatorial plane, 
the first one can be trivially chosen as 



si = dz = - - 



(4.32) 



The second one is obtained from the orthogonality and 
normalization constraints, and reads 



S2 — 

N 



A{r) 



dt + dr 



+ 



1 



&A(r) 

where the normalization function is 



(4.33) 



N 



bA{r) 



\T^l-A{r)Jl 



h^A{r) 



(4.34) 



Using the Sachs basis defined by Eqs. (4.30), (4.31), (4.32), 



and get 



and (4.33 ), we can finally compute the optical tidal matrix. 



-n{r) 

n{r] 



(4.35) 



where the function TZ{r) is 

3 /L 

As expected from the general decomposition (3.8), IZ is 



(4.36) 



trace free because only Weyl focusing is at work. Let us 
finally emphasize that A does not appear in the expres- 
sion (4.36) of 7?.(r), which is not surprising since a pure 



cosmological constant does not deflect light. 
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Wronski matrix 



where the functions B and C are given by 



The Sachs equations can now be integrated in order to 
determine the scale matrix Ck and the Jacobi matrix 1?k 
that compose the Wronski matrix Wk- 



First, since 7?. is diagonal, the Sachs equations (3.7) 



only consist of the following two decoupled ordinary dif 
ferential equations 



+n[riv)]Uv)- 



(4.37) 
(4.38) 



Clearly, the decoupling implies that the off-diagonal terms 
of Ck and f k vanish. 



^12 — ^21 — ^1^2 — ^21 — 0- 



(4.39) 



The calculation of the diagonal coefficients requires to 
integrate Eqs. (4.371 and (4.381. This cannot be per- 



formed analytically because there is no exact expression 
for r as a function v along the null geodesic. Indeed, we 
can write w as a function of r from Eq. (4.19) but this 
relation is not invertible by hand. 

Nevertheless, we are able to perform the integration 
perturbatively in the regime where £2/ei ^ ei <C 1, the 
relevance of which shall be justified by the orders of mag- 
nitude discussed in the next section. Solving Eq. (4.19) 
at leading order in £i, 62 leads to 



u{v) 



£1 



2 ■ (4-40) 



with Aw = h/E, and where denotes the value of the 
affine parameter v at the point of minimal approach. 



Eq. (4.37) then becomes, at leading order in ei,£2, and 



using the dimensionless variable w = (v — Wni)/Aw, 



dw'^ 



3£i / 1 
'2 ll + u'^ 



5/2 



6- 



(4.41) 



The perturbative resolution of Eq. (4.41) from v^^it to v 
finally leads to 



'--11 — ^ 



3ei 



- B'{winit){w - Winit) 



+ B{w) - B{w,^^) (4.42) 



and 



^ - "f^init) + ^ Aw |winit [B{w) ~ B{winit) 
- B'{winit){w - Winit)] - [C{w) - C(Winit) 

- B'{winit){w - Winit)] }+0 (^el , (4.43) 



1+9 7/)^ —7/1 

B{w) = ^^= and C{w) = —^=. (4.44) 



3VT 



The expressions of and I?^ are respectively obtained 
from Eqs. (4.42) and (4.431 by turning £1 into —ei. 

Note that in the limit £i,e2/£i — >■ 0, i.e. & — ?> 00 and 
A = 0; we find C = 1 and "D — (v — Winit)!, which are the 
expected expressions in Minkowski spacetime. 



C. Practical implementation 

This section has described the complete resolution of the 
equations for light propagation in a Swiss-cheese universe. 
All the results are included in a Mathematica program 
OneHole[] which takes, as input, the observation condi- 
tions and the properties of the hole; and returns Sm, fcin 

and W(Wsourcc;i'obs) = VVK(win;Wout) ' WfL (-^out ; Wout) ■ 

For simplicity, this program has been written assuming 
that the FL region has Euclidean spatial sections (K — 0). 

Iterating OneHole [] allows us to propagate a light 
signal back to an arbitrary emission event. Eventually, 
the redshift z is obtained by comparing the wave vector 
at emission and reception; and the luminosity distance is 
extracted from the block f (^source; '^obs) of the Wronski 
matrix W(wsourco; Wobs), according to 



(1 + z) ^ v^det I> ( ^source ; Vohs ) ■ 



(4.45) 



Note finally that, when iterating DneHole [] , we must 
also rotate the Sachs basis (31,82), to take into account 
that the plane of motion differs for two successive holes. 



V. EFFECT OF ONE HOLE 

Our method is first applied to a Swiss-cheese with a 
single hole. The purpose is to study the effects on the 
redshift and luminosity distance - for the light emitted 
by a standard candle - due to the presence of the hole. 



A. Numerical values of the parameters 

We assume that the clump at the center of the hole 
represents a gravitationally bound structure, such as a 
galaxy, or a cluster of galaxies. For a typical galaxy, 
we choose M ^ 10^^ Mq, and fo r a typical cluster M ~ 
IO^^Mq. By virtue of Eq. (2.19), the corresponding hole 
radii are respectively rh ~ 1 Mpc and rh ^ 200 Mpc. 



As discussed in the introduction, the typical width of 
a light beam associated with supernovae observations is 
^ AU. For comparison, the typical interstellar distance 
within a galaxy is ~ pc. Hence, such light beams are sen- 
sitive to the very fine structure of the universe, including 
the internal content of galaxies. In view of modelling the 
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relevant inhomogeneity scale with a Swiss-cheese universe, 
the above orders of magnitude suggest that the point 
masses inside the holes should represent stars, so that 
the natural choice should be M ^ Mq. Unfortunately, 
we cannot afford to deal with a so fine description, for 
numerical reasons. 

Nevertheless, the relevance of galactic holes in particu- 
lar can be justified by a signal/noise argument. A galaxy 
is indeed bright enough to flood the signal of a SN located 
behind it^. We conclude that the observation of a SN 
is possible only if its light does not cross any galaxy. In 
other words, galaxies can be considered opaque - though 
in a spurious sense - and their fine structure does not 
need to be modeled in first approximation. 

As discussed in Refs. [JHISS], the opacity argument 
provides a cutoff for the set of allowed impact parameters. 
Namely, we only consider light beams so that 



b > b„ 



''phys , 



(5.1) 



where Tphys is the physical size of the "opaque" object at 
the center of the hole. For a galaxy rphys ~ lOkpc, and 
for a cluster Tphys ~ 1 Mpc. 

The above discussions show that the various distance 
scales involved in our model are separated. The result- 
ing hierarchy is depicted on Fig. [7] and the typical or- 
ders of magnitude are summarized in Table |lj The lat- 
ter includes the expansion parameters £i = rg/b and 
£2 = Arg/3 ~ (fs/rHubbie)^- Their values justify a pos- 
teriori the perturbative expansion performed in § |IVB 6 
where we assumed that 62/61 <C £1 <C 1. In fact, it is 



even possible to show from Eq. (2.19) that £2 ^ £ 



Type 


rs (pc) 


rphys (kpc) 


rh (Mpc) 


£1 £2 


galaxy 


10"^ 


10 


1 


10"** - 10"® 10"^^ 


cluster 


100 


1000 


200 


10"'^ - 10"* lO"^'' 



TABLE I: Typical orders of magnitude for galaxy-like (M 
10"Mq) and cluster-Uke {M ~ lO^^Af©) Swiss-cheese holes. 



region is therefore characterized by the Einstein-de Sitter 
(EdS) cosmological parameters 



= 1, 



f^A = 0. 



(5.2) 



The effect of the cosmological constant will be studied in 



detail in Sec. VII The value of the Hubble parameter is 



fixed to Ho^ hx lOOkm/s/Mpc, with h = 0.72. 



B. Setup 

In order to study the corrections to the redshift z and 
luminosity distance Dj^, due to the presence of the hole, 
we consider the situation depicted on Fig. [8j 




FL 



Kottler 



'phys ^ '"s ^ 




'phys ''h 



''Hubble 



^(FL) 
r~ -^source 



FIG. 8; Setup for evaluating the effect of one hole on the 
redshift and luminosity distance. 



FIG. 7: Geometry and hierarchy of distances for a typical 
Swiss-cheese hole. 

In this section and the next one, we temporary set for 
simplicity the cosmological constant to zero. The FL 



^ The typical absolute magnitude of a galaxy is —22, while for the 
brightest SN it only reaches —19.3. 



Our method is the following. We first choose the 
mass M inside the hole and the redshift Zgourco of the 
source. We then fix the comoving distance between the 
observer and the center of the hole, in terms of the cos- 
mological (FL) redshift z^^^"^ of the latter. To finish, we 
choose a direction of observation, defined by the angle f3 
between the line of sight and the line connecting the 
observer to the center of the hole. 

Given those parameters, the light beam is propagated 
(in presence of the hole) until the redshift reaches Zgourcc- 
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We obtain the emission event ^source 
distance Dl 



We then compute ziom-co and D^^^ 



and the luminosity 
by 



considering a hght beam that propagates from ^source to 
the observer without the hole (bottom panel of Fig.l8|. 



Corrections to the redshift 



We only consider directions of observation such that the 
light beam crosses the hole. Thus, /^min < /? < ,5max where 
/3min and /3,„ax depend on the physical cutoff Tphys, the 
radius of the hole, and its distance to the observer 
Those dependences can be eliminated by plotting 6z as 
a function of (/3 - /3min)/(/5max - /3min) instead of ^, as 
displayed on Fig. \T0\ 



1. Numerical results 

The effect of the hole on the redshift is quantified by 

7 - z(FL) 



5z 



,(FL) 



(5.3) 



where we used the short notation z instead of Zsourcc- 
Fig. |9] shows the evolution of 5z with /3, for Zgourco — 0.05 
and various hole positions and masses. We have chosen 
M ^ IO^^Mq because the effect is more significant and 
displays less numerical artifacts than for M ^ lO^^Af0. 
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FIG. 9: Relative correction to the redshift z, due to the hole in 
the line of sight, as a function of the direction of observation /3, 
for a source at 2aourco = 0.05. Top panel: the mass of the hole 
is A4 = 10^^ Mq, and three positions between the source and 



the observer are tested, z^^^ / 



= 0.1 (blue, dot-dashed). 



0.5 (purple, dashed), and 0.9 (red, solid). Bottom panel: the 
hole is at z^^^ = 0.5 Zsomcc and three values for the mass are 
ii5 3 (-i^jyg^ dot-dashed), 2 (purple, dashed). 



tested, M/IO^^Mq 
and 1 (red, solid). 



E 4 











" ~ Z/, — 0. 1 ^ Z ^\iYce 






~ ~ ^- li ~ ^ ^ source 






Z/j — 0.9 source 











0.0 



0.2 0.4 0.6 

( /3-y3™, )/( ^„„-/3„,i„) 



0.8 



1.0 



o 
X 



M 


= 10 M 







- - M 


= 2xlO'^M 




o 


■ - - M 


= 3xlO'^M 




o 



0.0 



0.2 



0.4 0.6 

3m,„ ) / ( ySm„x " 



0.8 



1.0 



FIG. 10: Same as Fig. |9] but plotted in terms of the centered 
and normalized observation direction (/3 — /9min)/(/3max — /^min). 

As expected, 6z tends to zero when /? approaches /3niax 
(light ray tangent to the hole boundary). We notice that 
6z does not significantly depend on the distance between 
the observer and the hole. However, the effect clearly 
grows with the mass of the hole. 



2. 



estimation of the effect 



The correction in redshift due to hole can be understood 
as an integrated Sachs- Wolfe effect (see e.g. Chapter 7 
of Ref. 1 ). As the boundary of the hole grows with 
time (see Fig. [6]), the light signal undergoes a stronger 
gravitational potential at entrance than at exit. That 
induces a gravitational redshift (Szgi-av which adds to the 
cosmological one, and reads 



l + 6z„ 



'^out 



A(ro 



A(n, 



(5.4) 
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The order of magnitude of (5zgrav can be evaluated as 
follows. Let Sr = rout — Tin be the increase of the radius 
of the hole between entrance and exit. The expansion 
dynamics implies Sr ~ y/WiAt, where At — tout— ^in ^ Tim 
Tout is the time spent by the photon inside the hole. Using 
Eq. (5.4 1, we conclude that 



^^rav ^ • (5-5) 

For M — IO^^AIq (cluster- like hole), the numerical values 
given in Table □ yield (Szgrav,max ~ 10"*^. This order of 
magnitude is compatible with the full numerical integra- 
tion displayed on Figs. [9| and 

Such an analytical estimate enables to understand why 
dz increases with M, that is, with the size of the hole. 
Indeed, the bigger the hole, the longer the photon travel 
time so that the hole has more time to grow, and finally 
A{rout) - A{rin) is larger. 



D. Corrections to the luminosity distance 

The effect of the hole on the luminosity distance can 
be characterized in a similar way by 



6D^ 



(FL) 



D 



(FL) 



(5.6) 



The associated results, in the same conditions as in the 



previous paragraph, are displayed on Figs. 11 and|12[ 

We notice that is maximum if the hole lies halfway 
between the source and the observer, which is indeed 
expected since the lensing effects scale as 



Da (observer, lens) x D^^ilens, source) 
Da (observer, source) 



(5.7) 



which typically peaks for zicns ~ -Zsourco/S. The maximal 
amplitude of the correction is of order 10~^, for masses 
ranging from lO^^M© to 3 x IQ^^Mq. Just as for the 
redshift, the effect increases with the size of the hole. 

Note that 5Di^ can be related to the relative magnifi- 
cation /i, frequently used in the weak-lensing formalism, 
and defined by 
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(FL)^ 



Da 



l + z 

1 + z(FL) 
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(FL) 



Dl 



(5.8) 



Hence, if the correction on z is negligible compared to 
the one of Da, then the relation between SDi^ and /i is 

1 



6Dt 



1. 



(5.9) 
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FIG. 11: Relative correction to the luminosity distance Dl, 
due to the hole in the line of sight, as a function of the direction 
of observation /3, for a source at ^source = 0.05. Top panel: the 
mass of the hole is M = 10^^ M©, and three positions between 



source 



0.1 



the source and the observer are tested, z^^ / z. 
(blue, dot-dashed), 0.5 (purple, dashed), and 0.9 (red, solid). 

Bottom panel: the hole is at ^j^^^'"' = 0.5 ^source and three 
values for the mass are tested, M/10^^Mq = 3 (blue, dot- 
dashed), 2 (purple, dashed), and 1 (red, solid). 



and SDi^ ~ 100 Sz. The same study for M - 10"Mq 
leads to similar results with Sz ~ 10""'^'' — 10^^. Therefore, 
the effects of a single hole seem negligible. 



VI. EFFECT OF SEVERAL HOLES 

We now investigate a Swiss-cheese model containing 
many holes arranged on a regular lattice. Again, in all 
this section, the cosmological parameters characterizing 
the FL region are the EdS ones. 



E. Summary 

The presence of a single hole between the source and the 
observer induces both a correction in redshift and luminos- 
ity distance. For a hole with mass M ^ IO^^Mq, the rel- 
ative amplitude of those corrections are Sz ^ 10~^ — 10~^ 



A. Description of the arrangement of holes 

1. Smoothness parameter 

The smoothness of the distribution of matter within 
a Swiss-cheese can be quantified by a parameter / con- 



15 



o 
X 



I 



/' 






/ 

1 

_ 1 






. / 
1 


- - Zh=OA XZ source 




1 

;i 


^ /; ~ ^ Z source 
Z ll ~ X Z source 


N 

s. 

\ 

N 

\ 

\ 

\ 

\ 






T 




\ 



0.0 0.2 0.4 0.6 0.8 

(,6-,emm)/(,6m.-,x-jSmi„) 



1.0 



o 

X 



1 




M 


= 10'^ M 




o 


- - M 


= 2xlO'^M 







■ - - M 


= 3xlO'^M 








0.0 0.2 0.4 0.6 0.8 



1.0 



FIG. 12: Same as Fig. [TT] but plotted in terms of the centered 
and normalized observation direction (/? — /?min)/(Anax — /Smin)- 



structed as follows. Choose a region of space with - comov- 
ing or physical - volume V, where V^^^ is large compared 
to the typical distance between two holes. Thus, this 
volume contains many holes, the total volume of which 
is ViioiGs, while the region left with homogeneous mat- 
ter occupies a volume Vpi, = V — T^ioies- We define the 
smoothness parameter by 



lim — — . 

V-i-oo V 



(6.1) 



In particular, f — I corresponds to a Swiss-cheese with no 
hole ~ that is, perfectly smooth - while f — corresponds 
to the case where matter is under the form of clumps. Of 
course, / also characterizes the ratio between the energy 
density of the continuous matter and the mean energy 
density. 



of the smoothness parameter is in this case 

TT 



fm'in 1 



3\/2 



0.26. 



(6.2) 



In order to reach a smoothness parameter smaller than 
/mill, one would have to insert a second family of smaller 
holes. By iterating the process, one can in principle make 
/ as close as one wants to zero. 




FIG. 13: Hexagonal lattice of identical holes. On the left, the 
arrangement is close-packed, so that the smoothness parameter 
is / = /min « 0.26. On the right, / = 0.7. 



B. Observations in a unique line of sight 

We now focus on the corrections to the redshift and 
luminosity distance of a source whose light travels through 
the Swiss-cheese universe described previously. We study 
the influence of (a) the distance between the source and 
the observer, (b) the smoothness parameter /, and (c) 
the mass M of the holes. 



1 . Setup 

After having chosen the parameters (/, M) of the model, 
we arbitrarily choose the spatial position of the observer 
in the FL region, and fix its direction of observation. 
The method is then identical to the one of Sec. |Vl The 
light beam is propagated from the observer until the 
redshift reaches the one of the source, z. The ending 
point defines the emission event ^source- We emphasize 
that only emission events occurring in the FL region are 
considered in this article. 



2. Lattice 

We want to construct a Swiss-cheese for which the 
smoothness parameter is as small as possible. If all holes 
are identical, this close-packing problem can be solved by 
using, for instance, a hexagonal lattice. The corresponding 



arrangement is pictured on Fig. 13 The minimal value 



2. Influence of the smoothness parameter 

In this paragraph, the mass of every hole is fixed to 
M — IO^^Mq (galactic holes). The relative corrections to 
the redshift Sz and luminosity distance SDi^, as functions 
of the redshift z of the source, have been computed and are 
displayed on Fig. [14] for different values of the smoothness 
parameter /. 
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FIG. 14: Relative corrections to the redshift z (top panel) 
and luminosity distance Dl (bottom panel) as functions of z, 
for an arbitrary light beam travelling through a Swiss-cheese 
universe. All holes are identical, their mass is M — VS^^Mq. 
Three different smoothness parameters are tested: / = 0.26 
(blue, dot-dashed), 0.6 (purple, dashed), and 0.9 (red, solid). 



FIG. 15: Same as Fig. [Tijbut with / = 0.26 and three different 
values for the masses: M/10"Mq = 1 (blue, dot-dashed), 2 
(purple, dashed), and 3 (red, solid). 



Statistical study for random directions of 
observation 



While the corrections to the redshift remain small - 
typically bz < 10~^ - the cumulative effect of lensing 
on the luminosity distance is significant. For instance, a 
source at z 1.5 would appear 10% farther in a Swiss- 
cheese with / = 0.26, than in a strictly homogeneous 
universe. Both Sz and (5Z?l increase with z and decrease 
with /, as intuitively expected. Thus, the more holes, the 
stronger the effect. As examples, the light beam crosses 
~ 300 holes for (/ = 0.26, z = 0.1) or (/ = 0.9, z = 1), 
but it crosses ~ 2000 holes for (/ = /min, z = 1). 



3. Influence of the mass of the holes 

We now set the smoothness parameter to its minimal 
value /min ~ 0.26, and repeat the previous analysis for 
various hole masses. The results are displayed on Fig. 15 

We conclude that neither 6z nor i5Dl depend signifi- 
cantly on M, that is, on the size of the holes. Thus, what 
actually matters is not the number of holes intersected 
by the beam, but rather the total time spent inside holes. 



The previous study was restricted to a single line of 
sight, but since a Swiss-cheese universe is not strictly 
homogeneous, the corrections to z and Dl are expected 
to vary from one line of sight to another. As pointed out 
by e.g. Refs. [S31[57], such a restrictive analysis can lead 
to overestimate the mean corrections induced by inhomo- 
geneities. Besides, as stressed by Ref. [8], the dispersion 
of the data is crucial for interpreting SNe observations. 
Hence, the conclusions of the previous subsection need 
to be completed by a statistical study, with randomized 
directions of observation. 

Since the effect on the redshift is observationally negligi- 
ble, we focus on the luminosity distance. After having set 
the parameters (/, M) of the model, we fix the position of 
the observer in the FL region. Then, for a given redshift 
z, we consider a statistical sample of A'obs randomly dis- 
tributed directions of observation d E S^, and compute 
SDi,{z, d) for each one. 

Figure [T6| shows the probability distribution of SDi^ for 
sources at redshifts z = 0.1 (top panel) and z = 1 (bottom 
panel). We compare two Swiss-cheese models with the 
same smoothness parameter / = /min but with different 
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values for the masses of their holes {M = 10 Mq and 
IO^^Mq). The histograms of Fig. 16 are generated from 



statistical samples which contain A'obs = 200 directions 
of observation each. 
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FIG. 16: Probability distribution of the relative correction to 
the luminosity distance for random directions of observation, 
at z = 0.1 (top panel) and 2 = 1 (bottom panel). The 
smoothness parameter is / = /min and two different values for 
the masses of the holes are tested: M = lO^^Af0 (blue, solid) 
and M = lO^^M© (purple, dashed). 



From the statistical samples, we can compute the mean 
correction {6Di^)(z) and its standard deviation asD-^iz), 
whose evolutions are plotted on Fig. [iT) 

The results displayed on Fig. [TT] confirm the conclu- 
sions of Subsec. IVIBl The distance-redshift relation in 
a Swiss-cheese is biased with respect to the one of a 
purely homogenous universe. This effect is statistically 
significant, we indeed estimate (empirically) that 



(6.3) 



The bias slightly decreases with the mass parameter M. 
However, it can be considered quite robust because a 
variation of 4 orders of magnitude for AI only induces a 
variation of ~ 10% for the bias. 

The intrinsic dispersion of Z^l, associated with crsDi,, 
can be compared with the typical dispersion of the obser- 
vation. For instance, at z = 1 the former is ~ 1 %, while 




FIG. 17: Evolution, with redshift z, of the relative correction 
to the luminosity distance averaged over A'obs = 200 random 
directions of observation. Error bars indicate the dispersion 
o"i5Dl around the mean correction {SD-l). As in Fig. [Tm we 
compare Swiss-cheese models with M = IO^^Mq (blue, filled 
markers) and M — W^^Mq (purple, empty markers). 



the latter is estimated to be typically ~ 10% [67]. It 
follows that the dispersion induced by the inhomogeneity 
of the distribution of matter remains small compared to 
the observational dispersion. 



D. Summary and discussion 

This section has provided a complete study of the effect 
of inhomogeneities on the Hubble diagram, investigating 
both the corrections to the redshift and luminosity dis- 
tance of standard candles. The Swiss-cheese models are 
made of identical holes, defined by their mass M, and 
arranged on a regular hexagonal lattice. The fraction of 
matter remaining in FL regions defines the smoothness 
parameter /. For the hexagonal lattice, /min ~ 0.26. 

The effect on the redshift is negligible {Sz < 10~^), 
while the correction to the luminosity distance is signif- 
icant ((5Z3l > 10 % at high redshift). Compared to the 
homogeneous case, sources are systematically demagnified 
in a Swiss-cheese universe. The effect increases with z 
and decreases with /. 

Our results differ from those obtained in Swiss-cheese 
models with Lemaitre-Tolman-Bondi (LTB) solutions in- 
side the holes. In the latter case, a source can be either de- 
magnified if light mostly propagates through underdense 
regions [50l IHl [55] , or magnified otherwise. It has been 
proven in Refs. |S31 [571 [SB] that the global effect averages 
to zero when many sources are considered. Hence, LTB 
holes introduce an additional dispersion to the Hubble 
diagram, but no statistically significant bias. On the con- 
trary, in the present study, light only propagates through 
underdense regions, because we only consider light beams 
which keep far from the hole centers. This assumption 
has been justified in Subsec. [VA| by an "opacity" argu- 
ment. The bias displayed by our results is mostly due to 
the selection of the light beams which can be considered 
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observationally relevant. 

Our results also differ qualitatively from those obtained 
in the framework of the perturbation theory. In Ref. [55] . 
the probability density function of the weak lensing 
magnification fi, due to the large scale structure, has 
been analytically calculated by assuming an initial power 
spectrum with slope n — —2. Just as for LTB Swiss- 
cheese models, the magnification shows no intrinsic bias 
(i.e. (/i) = 1), but it is shown that P{fJ.) peaks at a 
value /ipcak slightly smaller than 1. Hence, a bias of 
order fipcak ~ (fj), which is typically 1% at z = 1, can 
emerge from observations because of insufficient statistics. 
However, this bias is far smaller than the one obtained in 
our Swiss-cheese model, of order 2SDi^ ^15%atz = l. 

Besides, the dispersion around the mean magnification 
is stronger for perturbation theory (~ 10 %) than for both 
LTB and Kottler Swiss-cheese models (~ 2%). 



VII. COSMOLOGICAL CONSEQUENCES 

Since the Hubble diagram is modified by the presence 
of inhomogeneities, the resulting determination of the 
cosmological parameters must be affected as well. 

More precisely, consider a Swiss-cheese universe whose 
FL regions are characterized by a set of cosmological pa- 
rameters (fiin, ^^a), called background parameters in 
the following. If an astronomer observes SNe in this inho- 
mogeneous universe and constructs the resulting Hubble 
diagram, but fits it with the usual FL luminosity-redshift 
relation - that is, assuming that he lives in a strictly ho- 
mogeneous universe -, then he will infer apparent cosmo- 
logical parameters which shall differ from 
the background ones. Evaluating this difference is the 



purpose of Subsecs. VHA VHB and VHC 



The natural question which comes after is: assuming 
that our own universe is well-described by a Swiss-cheese 
model, what are the background cosmological parameters 
which best reproduce the actual SNe observations? This 
issue is addressed in Subsec. IVII Dl 



A. Generating mock Hubble diagrams 

The Hubble diagram observed in a given Swiss-cheese 
universe is constructed the following way. We first choose 
the parameters of the model: /, M, and the background 
cosmology (Om, f^A = 1 — ^^m)'^- We then fix arbitrarily 
the position of the observer in the FL region, and we 
simulate observations by picking randomly the line of 



^ Recall that in the practical implementation of the theoretical 
results (see Subsec. |IV C i, we assumed that = 0, so that the 



background cosmology of our Swiss-cheese models is completely 
determined by f2m or SIa- Nevertheless, the apparent curvature 
parameter fix is a priori non zero. 



sight d, the redshift z € [0, Zmax], and we compute the 
associated luminosity distance Di^{z,d) as in Sec. VI In 
order to make our mock SNe catalog resemble the SNLS 3 
dataset [69], we choose Zmax = 1-4 and TVobs = 472. 

An example of mock Hubble diagram, corresponding 
to a Swiss-cheese model with / = /min, ^/ = 10^^ Mq 
and (fim, f^K, ^^a) = (1,0,0) is plotted on Fig. 18 As a 



comparison, we also displayed -Dl(^) for a homogeneous 
universe with (1) the same cosmological parameters, and 

(2) with {nrn,i^K,^A) = (0.3,0,0.7). 
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FIG. 18: Hubble diagram of a Swiss-cheese universe (dots) 
with / = /min, Af — IO^^Mq and EdS background cosmology. 
For comparison, we also display the distance-redshift relations 
of purely FL universes, with EdS parameters (blue, solid) and 
(nm,nA',nA) = (0.3,0,0.7) - (black, dashed). 



B. Determining apparent cosmological parameters 

The apparent cosmological parameters Hm, Ha and 
= 1 — ^^m — i^A are determined from the mock Hubble 
diagrams by performing a x^"fit- The is defined by 



472 r 



Ml — MFl(2i I ^^m, ^^a) 



(7.1) 



where [l no longer denotes the magnification, but rather 
the distance modulus associated with /^l, so that 



= 5 log 



10 



10 pc 



(7.2) 



In Eq. ( 7.1 1, (zi, /li) is the ith observation of the simulated 
catalog. In order to make the analysis more realistic, 
we have attributed to each data point an observational 
error bar A/x^ estimated by comparison with the SNLS 
3 dataset [BS]. Besides, /iFL(^ | ^m, ^^a) is the theoretical 
distance modulus of a source at redshift z, in a FL universe 
with cosmological parameters fim? ^Aj = 1 — Hni — f^A- 
The results of this analysis for two mock Hubble dia- 



grams are shown on Fig. 19 An EdS background leads 



19 



to apparent parameters (f^m, ^^x, ^^a) = (0.5,0.8,-0.3), 
which are very different from (1,0,0). Thus, the positive 
shift of Di^{z) - clearly displayed on Fig. [18]- turns out 
to be mostly associated to an apparent spatial curvature, 
rather than to an apparent cosmological constant. The ef- 
fect is weaker for a background with (f^m, ^^a) — (0.3, 0.7), 
which leads to {n,^,nK,^A) = (0.2,0.2,0.6). 
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FIG. 19: Comparison between background parameters 
(crosses) and apparent parameters (disks) for two Swiss- 
cheese models with / = fmin and M = W^^AIq. In blue, 
{n^,nA) = (1,0) leads to {n,n,^A)_= (0.5,-0.3). In black, 
(fin^^A) = (0.3,0.7) leads to (^^.^a) = (0.2,0.6). The la 
and 2a contours are respectively the solid and dashed lines. 
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FIG. 20: Apparent cosmic parameters f2m (gray disks), flK 
(red squares) and Hm (black diamonds) versus smoothness 
parameter /, for a Swiss-cheese universe with EdS background 
(Oni, flK, ^a) = (1, 0, 0). Solid lines and filled markers corre- 
spond to M = IO^^Mq, dashed lines and empty markers to 
M = lO^^M©. 
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C. Quantitative results 



Influence of the smoothness parameter 



FIG. 21: Apparent deceleration parameter g as a function 
of smoothness parameter /, for Swiss-cheese models with 
EdS background {q = 1/2). Solid lines and filled markers 
correspond to M = 10^^ Mq, dashed lines and empty markers 



Consider a Swiss-cheese model with EdS background 
cosmology. Fig. [20] shows the evolution of the apparent 
cosmological parameters with smoothness /. As expected, 
we recover fli = $7^ when / = 1, the discrepancy between 
background and apparent cosmological parameters being 
maximal when / = fmin- Surprisingly, a Swiss-cheese 
universe seems progressively dominated by a negative 
spatial curvature for small values of /. 

The apparent deceleration parameter q — f2,n/2 — Ha is 
plotted on Fig. 21 as a function of /. Interestingly, even for 
/ — /min, 9 remains almost equal to its background value 
q = 1/2. Therefore, though the apparent cosmological 
parameters can strongly differ from the background ones, 
the apparent expansion history is almost the same - at 
second order - as the background one. 



consistent with each other. The apparent cosmological 
constant f2A is slightly smaller for M = IO^^Mq than for 
M = 10^^ Mq, so that q is slightly larger. 



2. Influence of the background cosmological constant 

Now consider a Swiss-cheese model with / — /min and 
change its background cosmology. Fig. [22] shows the evo- 
lution of the apparent cosmological parameters versus the 
background cosmological constant f^A. As it could already 



Note that the results displayed on Figs. 20 and 21 are 



been suspected from Fig. [19] the difference between appar- 
ent and background parameters decreases with f^A, and 
vanishes in a de Sitter universe. This can be understood 
as follows. The construction of a Swiss-cheese universe 
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consists in changing the spatial distribution of the pres- 
sureless matter, while the cosmological constant remains 
purely homogeneous. Thus, the geometry of spacetime is 
less affected by the presence of inhomogeneities if fiA/f^m 
is greater. In the extreme case (fini,f^A) = (0,1), any 
Swiss-cheese is identical to its background, since there is 
no matter to be reorganized. 
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FIG. 22: Difference between apparent and background cos- 
mological parameters Clm — (gray disks), CIk — i^K (red 
squares) and Q,a^O,a (black diamonds) versus background Qa, 
for Swiss-cheese models with / = /min. Solid lines and filled 
markers correspond to M = IO^^Mq, dashed lines and empty 
markers to M = IO^^Mq. 



We can also plot the difference between the apparent 
deceleration parameter q and the background one q — 
f2m/2 — r^A, as a function of q. Again, q does not signifi- 
cantly differ from q. This result must be compared with 
Fig. 11 of Ref. gni, where {q-q)/q « 100%. 
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FIG. 23: Difference between apparent and background decel- 
eration parameters g — g as a function of q, for Swiss-cheese 
models with / = /min. Solid lines and filled markers corre- 
spond to M = 10^^ Mq, dashed lines and empty markers to 
M = IO^Mq. 



3. Comparison with other recent studies 

The impact of a modified luminosity-redshift relation 
~ due to inhomogeneities - on the cosmological param- 
eters had already been investigated by several authors. 
In Ref. [51], it has been suggested that a Swiss-cheese 
model with LTB holes and EdS background displays an 
apparent cosmological constant Ha = 0.4; but as already 
mentionned in Subsec. VI D[ su ch a claim was proven to 
be inaccurate in Refs. |53| 1541 [58] . because it relies on 
observations along a peculiar line of sight. When many 
random directions of observation are taken into account, 
the mean magnification goes back to 1. Hence, contrary 
to our results, the apparent cosmological parameters of a 
Swiss-cheese model with LTB holes are identical to the 
background ones. This conclusion is in agreement with 
Ref. |68j . where a similar study was performed in the 
framework of the perturbation theory. 

However, it is crucial to distinguish those approaches 
(LTB Swiss-cheese models and perturbation theory) to 
the one adopted in this article, because they do not ad- 
dress the same issue. The former share the purpose of 
evaluating the influence of inhomogeneities smoothed on 
large scales, while we focused on smaller scales for which 
matter cannot be considered smoothly distributed. Thus, 
our results must not be considered different, but rather 
complement ary. 



D. An alternative way to fit the Hubble diagram 

Let us now address the converse problem, and determine 
the background cosmological parameters of the Swiss- 
cheese model that best reproduces the actual observations. 
For that purpose, the simplest method would be to fit our 
observed Hubble diagram using the theoretical luminosity- 
redshift relation Dfp{z) of a Swiss-cheese universe. Hence, 
we need to derive such a relation in order to proceed. 



1. Analytical estimation of the distance-redshift relation of a 
Swiss-cheese universe 

As argued in Sec. |V) any observationally relevant light 
beam which crosses a Kottler region has an impact param- 
eter b much larger than the Schwarzschild radius rs of the 
central object. Moreover, since the cosmological constant 
has no effect on light focusing, we conclude that inside a 
hole, the evolution of the cross-sectional area of a light 
beam behaves essentially as in Minkowski spacetime. This 
conclusion is supported by the perturbative calculation 
of the Wronski matrix Wk performed in § |IV B 6| 

If both the observer and the source are located on 
the surfac e of a hole, their angular distance is therefore 
£)hoie _ y'(jet2> M Vout — I'm: wherc V denotes the affine 
parameter. More generally, for a beam which crosses N 
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contiguous holes, we get 



N 



7-^ holes 



5: A. 

i=l 



(7.3) 



where Avi = Wout,i — Vin,i is the variation of the afhne 
parameter between entrance into and exit from the ith 
hole. Let us now evaluate A?;^. The time part of the 
geodesic equation in Kottler geometry yields 



dt 



E r^rs 



dv A{r) 



E ~ constant, 



(7.4) 
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where E is the usual constant of motion. We co nclud e that 
« kl ^iAti- Besides, the relations (2.161 and (4.11) 



between FL and Kottler coordinates on the junction hy- 
persurface, together with ^(^h) ~ 1, lead to Ati « ATi 
flout /ao- Finally, 
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FIG. 24: Comparison between the approximate luminosity- 
redshift relation Dfp{z) = (1 + z)^D^{z) in a Swiss-cheese 
universe (solid lines), simulated observations (dots), and the 
Dyer-Roeder model D^^{z) with a{z) = / (dashed lines). 
Three different values of the smoothness parameter are tested, 
from top to bottom: / = fmin ~ 0.26, / = 0.7, / = 0.9. 



2. Comparison with the Dyer-Roeder approach 



where we approximated the sum over i by an integral. This 
operation is valid as far as ATi remains small compared 
to the Hubble time. In terms of redshifts, we have 



dz' 



(1 + z')2i7(z')' 



(7.6) 



By construction, this formula describes the behaviour 
of the angular distance when light only travels through 
Kottler regions. In order to take the FL regions into 
account, we write the distance-redshift relation D^{z) 
of the Swiss-cheese as the following (heuristic) linear com- 
bination 



Another widely used approximation to model the prop- 
agation of light in underdense regions has been proposed 
by Dyer and Roeder [T^ in 1972. It assumes that (1) the 
Sachs equation and the relation v{z) are the same as in 
a FL spacetime - in particular, the null shear vanishes 
- and (2) the optical parameter $00 



see 



IIIB2 



IS 



replaced by a(z)<I>00j where a{z) represents the fraction 
of matter intercepted by the geodesic bundle. In brief, 
the DR model encodes that light propagates mostly in 
underdense regions by reducing the Ricci focusing, while 
still neglecting the Weyl focusing. Under such conditions, 
the DR expression of the angular distance D^^{z) is 
determined by 



Dl^{z) = {l 



f) z?r°^(^)- 



(7.7) 



whe re / sti ll denotes the smoothness parameter defined 



VI A 1 and Dj^ (z) is the distance-redshift relation 



in a FL universe, given by Eq. (4.101 



A comparison between the above analytical estimation 
and the numerical results is plotted on Fig. [24j The agree- 
ment is qualitatively good, especially as it is obtained 
without any fitting procedure. Moreover, it is straight- 
forward to show that Dj^{z) and D^{z) are identical 
up to second order in z. This is in agreement with - and 
somehow explains - the numerical results of § VII C 1 and 



VII C 2 where we showed that the apparent deceleration 



parameter q is the same as the background one q. 

Note finally that the general tendency of our analytical 
relation is to overestimate (5Dl at high redshifts. The 
main reason is that its derivation uses the behaviour of 
Wk at zeroth order in rs/b, that is, it neglects the effect 
of the central mass in the Kottler region. The first-order 
term in Wk ~ taken into account in the numerical results 
- tends to lower the associated luminosity distance. 



d^-D™ 
dz2 



dlnH 
dz 

3 fiinO 



Ho 
H 



l + z 
2 



dp^ 
dz 



{l + z)a{z)D^^iz). (7 



This attempt to model the average effect of inhomo- 
geneities, while assuming that the universe is isotropic 
and homogeneous, has been widely questioned [7DH75] 
and recently demonstrated to be inconsistent [8] . 

Interestingly, our estimation Dj^(z) of the distance- 
redshift relation in a Swiss-cheese universe reads 



d^ 
dz^ 



dlnH 

3 f^mO 



l + z 

Ho 
H 



dz 



{l + z)fDl^{z). 



(7.9) 



which is similar to Eq. (7.8) with a{z) — /, except that 
the right-hand side is proportional to instead of D^. 
Nevertheless, it turns out that such a difference has only 
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a very weak impact, in the sense that 

\z)^D^^{z), i.e. Dl''{z)^D^^{z), (7.10) 



T^SC, 



if a{z) = f. This appears clearly on Fig. [24j where the 
dashed and solid lines are almost superimposed. In fact, 
it is not really surprising, since both approaches rely on 
the same assumptions: no backreaction, no Weyl focusing 
and an effective reduction of the Ricci focusing. 

Note however that this approach models the effect of 
the inhomogeneities on the mean value of the luminosity 
distance but does not address the dispersion of the data. 



3. Fitting real data with Dff^{z) 

The modified luminosity- redshift relation Dfp{z) de- 
rived in the previous paragraph can be used to fit the ob- 
served Hubble diagram. We apply the same method as 



described in Subsec. VII B except that now (1) the triplets 
{zi, ^i, A^i) are observations of the SNLS 3 catalog [Ml, 
and (2) fiFi,{z \ 0,^1,^1 a) is replaced by fj,sc{z \ ^m, f), 
where the background curvature VIk is fixed to (so that 
Oa = 1 — f^m)- Hence, we are looking for the smoothness 
parameter /, and the background cosmological parame- 
ters, of the spatially Euclidean Swiss-cheese model which 
best fits the actual SNe observations. 

The results of the fit are displayed on Fig. 25 First 
of all, we note that the confidence areas are very stretched 
horizontally, so that the smoothness parameter / cannot 
be reasonably contrained by the Hubble diagram. There 
are two reasons for this. On the one hand, we know 
from § |VII C 2| that / has only a weak influence on the 
luminosity-redshift relation in a universe dominated by 
the cosmological constant, which is the case here ^ 
0.7 - 0.8). On the other hand, since Dl'^{z) and Dl^{z) 
only differ by terms of order z^ and higher, one would need 
more high-redshift observations to descriminate them. 
However, all the current SNe catalogs - including the 
SNLS 3 dataset - contain mostly low-redshift SNe. 

Besides, Fig. 25 shows that fixing a given value of / 
changes the best-fit value of In the extreme case 
of a Swiss-cheese only made of clumps (/ = 0) we get 
= 0.3, while in the FL case (/ = 1) the best value is 
r^m — 0.24, in agreement with Ref. [69j. Such a discrep- 
ancy, of order 20%, is significant in the era of precision 
cosmology, where one aims at determining the cosmologi- 
cal parameters at the percent level. 

Let us finally emphasize that such a fit is only indicative, 
because it relies on an approximation of the luminosity- 
redshift relation in a Swiss-cheese universe. 



VIII. CONCLUSION 

In this article, we have investigated the effect of the 
distribution of matter on the Hubble diagram, and on 
the resulting inference of the cosmological parameters. 
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FIG. 25: Fit of the Hubble diagram constructed from the 
SNLS 3 dataset (BS], by using the luminosity-redshift relation 
Df^'~^ (z \Q,in, f) of a spatially Euclidean Swiss-cheese model. 
The colored areas indicate (from the darkest to the lightest) 
the la, 2a and 3a confidence levels. 



For that purpose, we have studied light propagation in 
Swiss-cheese models. This class of exact solutions of the 
Einstein field equations is indeed very suitable, because 
it allows for a strongly inhomogeneous distribution of 
matter, and this inhomogeneity does not backreact on 
the global cosmic expansion. The latter is entirely gov- 
erned by the background cosmological parameters fim, 
characterizing the FL regions of the model. The inho- 
mogeneities are clumps of mass M, while the fraction of 
remaining fluid matter is / - called smoothness parame- 
ter. The Swiss-cheese models are therefore defined by two 
"dynamical" parameters (f^mjOA), and two "structural" 
parameters (/, M). 

The laws of light propagation in a Swiss-cheese universe 
have been determined by solving the geodesic equation 
and the Sachs equation. For the latter, we have intro- 
duced a new technique - based on the Wronski matrix - in 
order to deal more easily with a patchwork of spacetimes. 
Our results, mostly analytical, have been included in a 
Mathematica program, and used to compute the impact 
of the Swiss-cheese holes on the redshift and on the lu- 
minosity distance. For a light beam which crosses many 
holes, we have shown that the effect on the redshift re- 
mains negligible, while the luminosity distance increases 
significantly with respect to the one observed in a FL 
universe (i5Z?l ~ 10% for sources at z ~ 1), inducing a 
bias in the Hubble diagram. 

The consequences of the bias on the inference of the 
cosmological parameters have been investigated by simu- 
lating Hubble diagrams for various Swiss-cheese models. 
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and by fitting them with the usual FL luminosity-redshift 
relation. In general, the resulting "apparent" cosmologi- 
cal parameters are very different from the "background" 
ones which govern the cosmic expansion, but in a way 
that leaves the deceleration parameter unchanged. More- 
over, the discrepancy between apparent and background 
cosmological parameters turns out to decrease with A, 
and is therefore small for a universe dominated by the 
cosmological constant. Finally, we have derived an ap- 
proximate luminosity-redshift relation for Swiss-cheese 
models, which is similar to the one obtained following 
the Dyer-Roeder approach. Using this relation to fit the 
Hubble diagram constructed from the SNLS 3 dataset, 
we have found that the smoothness parameter cannot 
be constrained by such observations. However, turning 
arbitrarily / = 1 into / = has an impact of order 20 % 
on the best-fit value of fim, which is significant in the era 
of precision cosmology. 

Of course, our model is oversimplifying for various rea- 
sons. First, it does not take into account neither the 
complex distribution of the large scale structures, nor the 
presence of diffuse matter on small scales - such as gas 
and possibly dark matter. Second, it does not take strong 
lensing effects into account, assuming that clumps are 
"opaque" . We can conjecture that this overestimates the ac- 
tual effect of the inhomogeneities. Nevertheless, it shows 
that their imprint on the Hubble diagram cannot be ne- 
glected, and should be modeled beyond the perturbation 
regime. Note finally that several extensions are allowed by 
our formalism. For instance, we could introduce different 
kinds of inhomogeneities, in order to construct fractal 
structures for which the smoothness parameter is arbi- 



trarily close to 0. Additionally to the Hubble diagrams, 
we could also generate the shear maps of Swiss-cheese 
models, and determine whether their combination allows 
for better constraints on the various parameters. 

This work explicitly raises the question of the meaning 
of the cosmological parameters, and of whether the values 
we measure under the hypothesis of a pure FL background 
represent their "true" or some "dressed" values. Similar 
ideas have actually been held in a similar context [ZH [75] . 
We claim that the simplest Swiss-cheese models are good 
models to address such questions ~ as well as the Ricci- 
Weyl problem and the fluid approximation - with their 
own use, between the perturbation theory and TV-body 
simulations. 
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